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Abstract. We consider the XXX spin-^ Heisenberg chain on the circle with an arbi¬ 
trary twist. We characterize its spectral problem using the modified algebraic Bethe anstaz 
and study the scalar product between the Bethe vector and its dual. We obtain modified 
Slavnov and Gaudin-Korepin formulas for the model. Thus we provide a first example of 
such formulas for quantum integrable models without U(l) symmetry characterized by an 
inhomogenous Baxter T-Q equation. 
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The study of quantum integrable models with U(l) symmetry by the Bethe ansatz (BA) 
methods [4, 11, 34] provides exact solutions which found applications in a wide range of do¬ 
mains such as: statistical physics, condensed matter physics, high energy physics, mathematical 
physics, etc. One of the major accomplishments of the method has been the obtaining of form 
factors, for models related to 0(2 and 0(3 families of symmetry, in the compact form of deter¬ 
minants [8, 9, 21, 22, 24, 26, 31]. In particular, for models related to the 0(2 symmetry, the 
key results are the Slavnov [35] and the Gandin-Korepin [19, 20, 25] formulas, which provide, 
respectively, the scalar product between an eigenstate and an arbitrary state and the norm of 
the eigenstates. 

In the case of models without U(l) symmetry, the usual BA techniques in general fail to 
provide a complete description of the spectrum^. Thus alternative methods have been de¬ 
veloped, for instance, the separation of variables (SoV) [15, 16, 23, 30, 33], the commuting 
transfer matrices method [3], the functional method [18] or the g-Onsager approach [2]. Recently, 
key steps have been accomplished for the Bethe ansatz solution of such models. On the one 
hand, a new family of inhomogeneous Baxter T-Q equation to determine the eigenvalues has 
been proposed by the off-diagonal Bethe ansatz (ODBA) [12, 37]. On the other hand, the 
construction of the off-shell Bethe vector has been done in the context of a modified algebraic 
Bethe ansatz (MABA) approach [1, 5, 6, 10]. Let us remember that previous developments in 
the BA technique, in particular the obtainment of the eigenvectors of the XXX chain on the 
segment with upper-triangular boundaries [7], brought important insights to the MABA. 

^In some cases, some gauge transformation can allow to apply the ABA, see for example the XYZ spin 
chain [4, 36]. For the XXX case, that we consider here, the GL(2) symmetry allows one to restore the U(l) 
symmetry [14] and the usual ABA applies, provided that the twist is a non-singular matrix. Also, in the context of 
open XXZ spin chains, constraints on the parameters of the model allow one to apply the usual techniques, see [28]. 
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Here, we consider the question of the scalar product between the Bethe vectors obtained from 
the modified algebraic Bethe ansatz. For simplicity, we consider the case of the isotropic spin—^ 
Heisenberg chain on the circle with an arbitrary twist. This is the simplest model which can be 
considered by the MABA. In fact, its solution contains the main features entailed by the method: 
the spectrum is characterized by an inhomogeneous Baxter T-Q relation and the off-shell Bethe 
vector is generated by a modified creation operator. In this context, we obtain a modified 
Slavnov formula for the scalar product between an on-shell Bethe state and its off-shell dual. 
The formula (see (17)) is given in terms of a determinant depending on the Jacobian of the 
inhomogenous eigenvalue expression. Moreover, it contains a new factor related to a certain 
expansion of the Bethe vector. The square of the norm, i.e., the modified Gaudin-Korepin 
formula (see (19)), is obtained by a limit. As expected, the case with general integrable open 
boundary, which also breaks the U(l) symmetry, has the same structure and will be presented 
in a separated publication. 

The isotropic spin—^ Heisenberg chain on the circle with an arbitrary twist is given by the 
Hamiltonian 

N 

^ C^k+1 +^l® ^l+l +<^l® (1) 

k=l 


subject to the following boundary conditions 


l^N+l = 


- k\- .K? - y _ 

-F Z ————- - Uf -k [KK- - KK+)af, 


y + K_ R +K +hi+ + n_ y 

TO-^+i = *■—-^^- -(^1 - 

= (KK+ — KK-)af + i{KK- + KK+)af + {kK + K+K_)(Tf. 


The twist parameters {k,k,k+,k-} G are generic and 7 = kk — k+k_. The Pauli matri¬ 
ces^ with a = x,y,z act non trivially on the fcth space of the quantum space 77 = 
with Pfc = C^. 

The Hamiltonian (1) is integrable and can be considered within the quantum inverse scatte¬ 
ring method [17]. Let us briefly recall this formalism. The key object is the rational i?-matrix 


R{u) = - + P, 
c 

which acts on ( 8 > C^, with P = ® a~ + a~ ® ® c G C*. From the 

7?-matrix, we construct the monodromy matrix 


which acts on ( 8 > 77 and with {0} = {0i,... 6j^'\ G C'^ being the inhomogeneity parameters. 
This monodromy matrix satisfies the RTT relation 



Rah{u - v)Ta{u)Ti,{v) = Tb{v)Ta{u)Rab{u - v) 

that encodes commutation relations between the operators {tij{u)}, see Appendix A. The trans¬ 
fer matrix, generating function of the conserved quantities of the model, is given by 


t{u) = Tra {KaTa{u)) = Ktii{u) -k Kt22{u) -k K+72 i(m) "k K 


0 -1 


T+ — 


0 1 
0 0 


0 0 
1 07 ’ 


= (T+ + a~ 


= i{a — I 
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with 


K = 


K 

K~ K 


( 2 ) 


The commutation relation between transfer matrices with different spectral parameters 
[t{u),t{v)] = 0, 

follows from the RTT relation and the GL(2) invariance of the i?-matrix 
[Rab{u - v),KaKb] = 0 . 

The Hamiltonian (1) is given by 


( 3 ) 


d 




'LL —>^0, Oi —^0 


-iV, 


and thus its spectral problem is the same of the one of the transfer matrix. 

The diagonalization of the transfer matrix can be obtained by means of the MABA [1, 5, 6, 10] 
and leads to an inhomogeneous Baxter T-Q equation. In order to do that, we introduce the 
following transformation of the twist matrix (2) 


K = LDL 


with 


r 1 M ^ 


D = 


R — p 0 
0 K — p 


= 


h + K — p 
R + K — 2p 


and 


p^ — [k + k)p + K = 0 . 

By means of this transformation, we can obtain a modified monodromy matrix T{u) = LT{u)L 
with entries given by modified operators {t'jj(u)}. They are expressed in terms of the initial 
{tij{u)} operators by 

( 2 

^ t22{u] 

K K~'~K 

l^uiu) = P (^uiu) + + t22iu)) + t2l{u)^ , 

= p (t2l{u) + + t22{u)) + tuiu)^ , 

1'22{u) = p (t22(u) + + -^t2l{u) + ^ _ tn{u)] ■ 

It follows that the transfer matrix has a modified diagonal form given by 
t{u) = Tr [Df{u)) = {k- p)un{u) + {k - p)i' 22 iu). 


To construct the Bethe vector we use the usual highest weight representation and the highest 
weight vector 
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The actions of the operators {tij{u)} on it are given by 
tii{u)\ 0 ) = Ai(tt)|0), t 2 i{u)\ 0 ) = 0, 

with 

^ a , ^ a 

i=l i=l 

For the dual Bethe vector we use the dual highest weight vector 
( 0 | = ( 1 , 0 )®", 
with the actions 


( 4 ) 


(0|tii(u) = Ai(u)(0|, (0|ti2(u) = 0, (0|0) = 1. 

For the operators {i'ij(u)}, we have a modified action on the highest weight vector. We can 
show that 


i^ii(u)|0) = Ai(rt)|0) + ■^i'i2{'u)\0), 
i^22iu)\0) = A2(it)|0) + -^J^i2{u)\0), 
J^2iiu)\0) = -^(Ai(tt) + A2 (u))|0) + 


1^12 (u) 10 ). 


( 5 ) 

( 6 ) 

( 7 ) 


We have thus all the ingredients to implement the MABA. We will use the notation u with 
= M for the set of M variables {ui,U2, ■ ■ ■ ,um}- If the element Ui is removed, we note 
Ui = {ui, U 2 ,.. ■, Ui-i,Ui+i ,..., um}- If we also remove the element Uj, we note Uij = u/{ui,Uj}. 
For products of functions (see (8) bellow) or of operators {i'ij{u)}, we use the convention 


M 


M M 

g{v,u) = n 

i=lj=l 
M 

g{ui,Ui)= g{ui,Uj), i^ijiu) = 


g{u,u) = Y\g{u,ui), 
i=l 

M 


k=l 


The functions 


g{u,v) = 


u — V 


f{u,v) = l + g{u,v) = 


u — V + c 


u — V 


( 8 ) 


will be widely used. 

Let us consider the vector 


sMi'- 


u) = Ui2{ui) ■ ■ ■ Ui2{um)\0) = Z/I 2 (n)| 0 ), 


(9) 


and act with the transfer matrix on it. In order to perform this calculation we need to find the 
action of the operators {r'jj(tt)} on the vector (9). Using the GL(2) invariance (3), it is easy 
to see that the new operators satisfy the same commutation relations (23), (24), (25) of the 
operators {tij{u)}, see Appendix A. Thus, using these commutation relations and the action 
(5), (6), (7), we can show that 
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{u) = + Xi{u)f{u,u)M^ (u) 

M 

+ ^ 9{u, Ui)\i{ui)f{ui, Uj), 

i=l 

i'22{u)'^^ {u) = + A2(n)/(u,'u)]B^(u) 

M 

+ ^ g{ui,u)X2{ui)f{ui, Ui)M^{u, Ui), 
i=l 

i/2i{u)B^{u) = u) 

/ M 


2 = 1 


M 


+ E^( 


sM-1 


u,Ui,Ui)Sii- ^(Ui)+ ^ G{u,Ui,Uj, 

i=l 


u. 


iM-l 


IJJ 


{u, u. 


( 10 ) 


• 13 / 5 


with 


A^(u, m|x, ?/) = x/(u, ?x)Ai(u) + yf{u, u)X2{u), 

E^{ui,Ui\x,y) = -xf{ui,Ui)Xi{ui) + yf{ui,Ui)X2iui), 

F{u, Ui, Ui) = g{u, Uj)Ai(?x)A2(ui)/(u, Ui)f{ui,Ui) + g{ui, u)Xi{ui)X2{u)f{ui, Ui)f{ui, u), 
G{u, Ui, Uj,Uij) = g{u, Ui)g{uj,u)Xi{ui)X2{uj)f{ui, Uj)f{ui, Uij) f (uij, uj) 

+ g{u, Uj)g{ui,u)Xi{uj)X2{ui)f{uj,Ui)f{uj,Uij)f{uij,Ui), 

and where we have used the functional identities 

M M 

f{u, u) + '^ g{u, Ui)f{ui, Ui) = I, f{u, F) + ^ g{ui, u)f{ui,Ui) = 1 . 

2=1 2=1 

It follows that the action of the transfer matrix is given by 

t{u)M^{u) = —M^~^^{u, u) + A^(tt, u\k — p,K, — p)E^{u) 

g 

M 

+ ^ g{ui, u)E^{ui,Ui\k - p,K- p)B^(m, Ui), ( 11 ) 

2=1 

where we have used the relation -^{h + K — 2p) = The new term B'^'’'^(n, u) is characteristic 
for models which break the U(l) symmetry. From the formula (11) we can obtain, by limit the 
upper triangular case = 0 and the diagonal case k~ = = 0. The Bethe ansatz solution is 

then obtained by requiring that E^{ui,Ui\k — p,K — p) = 0 for M = 0,..., N. Only the part 
with positive total spin contributes to the solution. Here to complete the MABA we must thus 
find the action of the operator 1^12 (w) on the Bethe vector (9) when M = N. In this case the 2^ 
independent diagonal Bethe vectors with M < N and partitions of N variables are contained 
in the vector, see Appendix B for the expression of the Bethe vector in terms of the initial 
operator ti 2 {u). This allows us to find an off-shell action with a wanted/unwanted form 

N 

—B^+^(u, n) = Ag {u,u)'E^{u) -\-'^g{ui,u)E^{ui,Ui)'E^{u,Ui), 
h- ,_i 


(12) 
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Ag {u, u) = 2p\i{u)\2{u)g{u, u), Eg {m, Ui) = 2pXi{ui)X2{ui)g{ui, Ui). 

This action can be proved following the method given in [13] and will be considered elsewhere. 
Finally, from (11) and (12), we obtain the off-shell equation with generic u with ffu = N, 

N 

t{u) M^{u) = A^{u, u)M^{u) + g{ui,u)E^{ui, Ui)M^{u, Ui), 

i=l 

with an inhomogeneous eigenvalue 

A^(m, u) = {R- p)Xi{u)f{u, u) + {k- p)X 2 iu)f{u, u) + 2pXi{u)X2{u)g{u, u), (13) 

and an inhomogeneous Bethe equation 

E^{ui,Ui) = -{R - p)Xi{ui)f{ui,Ui) + {k- p)X2{ui)f{ui,Ui) 

+ 2pXi{ui)X2{ui)g{ui,Ui), (14) 


for i = 1 ,..., A^. 

One can proceed in a similar way for the dual Bethe vector 

C^(u) = (0|i/2i(h), (15) 

and, in particular, obtain 

N 

C^{u)t{u) = {u,u)C^(u) + gjuj, u)E^{ui,Ui)C^{u, Uj), 

i=l 

with the same eigenvalue and Bethe equations of the Bethe vector (9) with M = N. 

When the Bethe equations are satisfied, i.e., E^{ui, Ui) = 0 for i = 1,N, and we consider 
non-singular solutions of the Bethe equations [29], the on-shell Bethe vectors are eigenstates of 
the transfer matrix 

[u] = A^ {u,u)'M^ {u), {u)t{u) = A^ {u,u)C^ (u). 

Let us remark that the completeness of the solution given by (13), (14) has been numerically 
checked for chains with small size. It should be interesting to prove it along the lines of [27]. 

We are now in position to consider the scalar products for the Bethe vector (9) and (15), 
namely, 

S^{u\v) = C^{u)M^{v). (16) 


From the construction given hereafter, we obtain the modified Slavnov formula of the twisted 
XXX spin chain characterized by the inhomogeneous Baxter T-Q equation (13). When 

E^{ui,Ui) = 0 , 


for i = 1,..., N, the scalar product (16) has a compact form given by 

,2 


S^{u, v) = c 


.N 


A 


R-\- K — p 


^ DetAT (J^A'^(uj,u)) 


DetAT [g{vi,Uj)) 


(17) 


with 


-\ N 


Wi\u)={ — ] {0Mu)\0), 
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which is given explicitly in Appendix B. When 

E^{vi,Vi) = 0 , 

for i = 1 ,..., A^, the scalar product (16) is given by 


S^{u, v) = c 




R + K — p 


Detiv {g{ui,Vj)) 


(18) 


We can obtain the square of the norm by imposing the limit u = v. Using the well-known 
formula for the Cauchy determinant 


Detiv {g{vi,Uj)) = 
we obtain 


g{v,u) 


n 9iui,Uj)g{vj,Vi)' 
i<j 


N 


N 


<{®) I n g{ui,uj)g{uj,ui) j DetAr(G- 

i<j 






+ K — p 

where the matrix elements Gij, for z, j = 1 ,..., A", are given by 

Gii = 2pc{\2{ui)dui\i{ui) + \i{ui)dui\2{ui)) 

( ^ 

+ {-l)^{R-p)\ch{ui,Ui)dui\i{ui)-\i{ui) ^ h{uij,Ui) 

( ^ 

+ {k- p) I ch{Ui,Ui)duiX2{Ui) + \2{Ui) 

Gij = {-1)^{R - p)Xi{uj)h{uij,Uj) - (k- p)X 2 {uj)h{uj,Uij) for i / j, 


(19) 


with 


h{u,v) = 


fiu,v) 

g(u,v} 


u — V + c 
c 


and where we have applied the I’Hospital’s rule to find 


Gij = lim c- 

Vj^Uj 


givj,u) 


( 20 ) 


From the conjectured modified Slavnov and Gaudin-Korepin formulas we remark that the 
ratio, needed for the calculation of the correlations functions, is independent of the Wq (m) and 
of the constant (Thus the relevant part of the formula for the correlations functions 
is given by the Jacobian of the inhomogeneous eigenvalue (13) and it limits (20). 

To obtain this conjecture we have proceeded in the following way. We start from the following 
hypothesis: 

• We can impose the Bethe equations (14), E^{vi, Vi) = 0 for i = 1,..., A^, by linearizing 
the quadratic terms Xi{vi)X 2 {vi) in terms of Ai(uj) and X 2 {vi). 

• The usual Slavnov formula must be recovered in the U(l) symmetric limit and thus the 
modified Slavnov formula contain the determinant of the Jacobian of the inhomogenous 
eigenvalue (13). 
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For the U(l) symmetric case {k~^ = k~ = 0), the scalar product between an on-shell Bethe 
vector and an off-shell Bethe vector is given by the Slavnov formula [35]. If the v are on-shell, 
i.e., E^{vi, Vi\K, k) = 0 for i = 1,..., M, we have 


S^{u,v) 


k) T)etM{g{ui.,Vj)) 


From the first term of action the (10), for M = N and using (12), we can extract the leading 
term of the scalar product with N quadratic terms: Xi{ui)X 2 {ui) with i = If we 

impose the Bethe equations (14), E^{vi,Vi) = 0 for i = 1,..., A", by linearizing the quadratic 
terms Ai(ui)A 2 (uj) in terms of Ai(ui) and X 2 {vi), this leading term is invariant and is the only 
contribution at the top order 3N in the A*. This term is proportional to 


N 

i=l 


It allows us to fix the functional to complete the determinant part of the Slavnov formula (18). 
Indeed, the functional is the same, up to a constant, when we consider the leading coefficient 
in the Aj and use the identity 


In the U(l) symmetry limit, p = 0, we have Wq {v) = A 2 (u)(f + 1)'^ and restore the usual 
Slavnov formula up to a constant. The formula (17) can de derived, also up to a constant, in 
a similar way starting from the action of the operator i>i 2 {u) on the dual Bethe vector (15). This 
way to find the modified Slavnov formula does not fix the constant. 

Let us fix this constant by considering the simplest case N = 1 and discuss the difficulties 
to find the good parametrization for the off-shell scalar product, which allows us to find the 
modified Slavnov formula by linearizing the quadratic term of the Bethe equation, for general N. 
In the case N = 1, the good one is 


r.N 


DetAT [g{uj,Vi)) {ui,v). 


2 = 1 



5^(u,u) = g f5](u,u) + ^:r^ -(A^(u,'(;)lTo^(u) -b A^(u,m)WoH^^)) 

V K, ~\~ Kj — p 


( 21 ) 


with 


Wq{u) = Xi{u) -b A 2 (u), S\{u,v) = g{u,v){Xi{v)X 2 {u) - Ai(u)A 2 (u)). 

To obtain the Slavnov formula we use the Bethe equation to linearize the quadratic term that 
depends on v 

Xi{v)X2{v) = ^((it - p)Ai(u) - (k - p)A2(u)) (22) 

and thus we arrive to our result 

^ ff(v,u)Wo\v)((K - p)Xi(u) - (k - p)X2(u) - 2pXi(u)X2(u)). 

If we also linearize the second quadratic term, we find zero for u ^ v, which shows the orthogo¬ 
nality of the Bethe vectors and the modified Gaudin-Korepin formula (19) for u = v. 
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The parametrization (21) is not unique; for instance, from the projection on the {tij{u)} 
operator, see Appendix B, we find 

u) = v) + ^ _ (Ai(u) + A2 (u)) (Ai(u) + Ai(u)) 

that is equivalent to (21) when we specify the explicit form of the Xi{u) (4). In this case the 
quadratic term Ai(u)A2(u) does not appear and the formula (22) could not be used directly. 
Moreover, the recursion relation (10) gives also another parametrization in order 3 in the Xi{u) 
like (21) but with terms of the form Xi{u)Xi{v). In this case the linearization of the quadratic 
term (22) could be used but does not lead directly to the modified Slavnov formula. A systematic 
way to fix the good parametrization of the off-shell scalar product (i.e., a parametrization that 
reduces to the modified Slavnov formula, when we linearize the quadratic terms Ai(ui)A2(uj) 
from the inhomogenous Bethe equations) remains an open problem for the moment and will be 
discussed elsewhere. The conjecture can be tested exactly for = 1, which allows us to fix the 
constant, and then the cases N = 2 and N = 3 have been checked numerically to support the 
conjecture. 

Finally, let us point out another problem which is worth to be considered: the construction 
of an explicit algebraic link between the MABA solution and the solution obtained by means of 
the usual ABA, see, e.g., [32]. One first step to address this problem, at least at the spectral 
level, is to equate the eigenvalue expression found from the MABA (13) with the ones obtained 
from the usual ABA. The simplest eigenvalue from the usual ABA is given by 

q;Ai(u) + {k + k — a)X2{u) 

with a = ^{k + k + {n — k)^ -|- 4 «:+k“) an eigenvalue of the twist matrix K. As a result, we 
can obtain constraints on the parameters u that, up to symmetrization, provide one solution of 
the Bethe equations (14). More details will be given elsewhere. 

A Rational functions and commutations relations 

Let us introduce the commutation relations of the Lj(u) given by 
[tij{u),tki{v)] = g{u,v)[tkj{v)tii{u) - tkj{u)tii{v)). 

In particular, we will only use the following ones 

= f{v, u)ti 2 {v)tii{u) + g{u, v)ti 2 {u)tii{v), (23) 

t 22 {u)ti 2 iv) = f{u, v)ti 2 {v)t 22 {u) + g{v, u)ti2(u)t22(t'), (24) 

t 2 liu)tl 2 {v) = ti 2 {u)t 2 l{v) + g{u,v)(tii {v)t22{u) - tu{u)t22{v)) . (25) 

The functions / and g are given by (8). They allow us to find the action on the multiple product 

of ti 2 {u) that forms a string of length M, namely, 

M 

tll{u)tl2{v) = f{v,u)tl2{v)tu{u) + g{u, Vi)f{vi, Vi)tl2{u)tl2{Vi)tll{Vi), 

i=l 

M 

t 22 iu)ti 2 {v) = fiu,v)ti 2 iv)t 22 {u) + g{Vi, u) f (Vi, Vi)ti 2 {u)ti 2 {Vi)t 22 iVi), 

i=l 

t2l{u)ti2{v) = ti2{v)t2l{u) 

M 

+ '^fiUi,Ui){g{u,Ui)t 22 {Ui)ti 2 iUi)tii{u) + gixi,y)t 22 {u)ti 2 {ui)tii{ui)). 
i=l 
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B Projection of the Bethe vector 

The Bethe vector in terms of the operator ti 2 {u) is given by 

N 

iyi2{u)\0) = ^^(«i|^n)tl2(Wn)|0), 

i=0 u^{ui,u]i} 

with = i, ^Ui = N — i a partition of the set u. The sum is over all ordered partitions, 
denoted u —>■ {'Ui,iin}. The coefficient is given by 

,... ,UAr)^ , 

where 


( N-i 


lT'i_i(ni|u2, ...,Ui) = f{ui,ui)Xi{ui) + /(tti,ui)A2('Ui) = ^{ui,ui, 1,1) 


and 

Symr(F(u)) = ^ F(r), 

o'^Sm 

with = {?Xo-(i), • • • ,Ufj(^M)} ™ element of the permutation group Sm- The dual Bethe vector 
in terms of the operator t 2 i{u) is given by 

N 

{0\ui2{u) = (^) Wl^{Uj\Uj,){0\t2l{Uu)- 

i=0 
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